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We give a survey on the predictions of Coulombic phase transitions in dense plasmas (PPT) and 
derive several new results on the properties of these transitions. In particular we discuss several 
types of the critical point and the spinodal curves of quantum Coulombic systems. We construct 
a simple theoretical model which shows (in dependence on the parameter values) either one alkali- 
type transition (Coulombic and van der Waals forces determine the critical point) or one Coulombic 
transition and another van der Waals transition. We investigate the conditions to find separate 
Van der Waals and Coulomb transitions in one system (typical for hydrogen and noble gas-type 
plasmas). The separated Coulombic transitions which are strongly influenced by quantum effects 
are the hypothetical PPT, they are in full analogy to the known Coulombic transitions in classical 
ionic systems. Finally we give a discussion of several numerical and experimental results referring 
to the PPT in high pressure plasmas. 



I. FIRST ORDER VAN DER WAALS AND 
COULOMB TRANSITIONS 

The theory of gases developed in the dissertation of 
van der Waals in 1873 may be considered as the starting 
point of the modern theory of phase transitions. Van der 
Waals approach is based on a simple physical model of 
interactions between particles which takes into account 
short range repulsive as well as long range attractive 
forces. The van der Waals model predicts below the 
critical temperature the possibility of coexistence of two 
phases which differ from each other by the density of 
molecules. In connection with the development of more 
strict theories it became clear that van der Waals ap- 
proach is restricted to relatively weak attractive forces 
which either decay faster than 1/r 3 or fulfil the Kac con- 
ditions. Therefore the applicability to Coulomb forces 
remained open. 

In 1943 Landau and Zeldovich discussed new possibil- 
ities of phase transitions connected with metal-insulator 
transitions in metal fluids in the vicinity of critical points 
jl| . We have to mention in this connection also the work 
of Mott on metal-insulator transitions in Coulomb sys- 
tems. A first systematic study of Coulombic transitions 
in plasmas was given by Norman and Starostin in 1968 
Q . The first order transition predicted by these authors 
was named plasma phase transition (PPT). The PPT was 
predicted as a possible result from the competition be- 
tween effective Coulomb attraction and quantum repul- 
sion in the partially ionized dense plasma |^|| . The qual- 
itative picture was similar to the van der Waals model 
where the phase transition is a result from the long-range 
attraction between neutral molecules and their short- 



range repulsion. We remember that the two phases in 
the van der Waals model differ from each other by the 
density of molecules. The two phases in PPT have dif- 
ferent number densities of charged particles and different 
degrees of ionization. Atoms, which are present in both 
coexistent phases, were treated originally as an ideal gas 

0- 

To get the first estimate of the critical temperature 
Norman and Starostin Q used the thermodynamic func- 
tions available at that time, namely Debye-Huckel ex- 
pressions for the chemical potential of the charges with 
quantum corrections The expressions described an 
effective Coulomb attraction and an effective quantum 
repulsion due to the uncertainty principle. In the linear 
approximation in the density a critical temperature T cr = 
2660 and with the more realistic nonlinear Debye-Huckel 
expression the value T cr = 10640 was obtained. We men- 
tion that both expressions are valid only in the limit of 
small quantum repulsion (X/rrj) << 1 (here A is the elec- 
tron thermal wavelength and ru is the Debye radius), 
and small non-ideality parameter 7 = e 2 ?!^ 3 / 'kT « 1, 
where rii is the ion number density, T is the temperature. 
Later estimates of the PPT in hydrogen plasmas lead to 
higher values of the critical temperature [5-14] 

The hypothetical phase transitions in multiple ionized 
plasmas were treated first in |lj|. A detailed study of 
He-plasmas was given in Jl6[ . Here we will consider only 
single-charged ions in gaseous plasmas. 

Since up to now the PPT in plasmas is not yet clearly 
identified experimentally we plan to make as clear as pos- 
sible the definition, the properties and the conditions for 
a PPT. For this reason we develop a simple model (a com- 
bination of the Debye-Huckel with the Van der Waals the- 
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ory) which shows a separate Coulombic transition - which 
is a PPT. Other related phenomena in plasma systems 
are discussed only in brief. 

Plasma-like phase transitions in the optically excited 
electron- hole system in semiconductors were studied with 
similar methods |^,^,|l^,^8|,|9| . The experimental obser- 
vations seem to point out that these instabilities of the 
theory correspond to a real phase transition in semicon- 
ductors, this topic will not be considered here in detail. 

In more detail we will discuss comparison of the PPT 
to the known transitions in classical Coulombic systems. 
The first results on phase transitions in ionic systems 
go back to the late sixties. In 1970 Voronzov et al. ]2C| ] 
found a coexistence line and a critical point in the course 
of Monte Carlo studies of charged hard spheres imbed- 
ded into a dielectric medium. An analytical estimate of 
the critical point and the coexistence line for electrolytes 
based on the Debye-Hiickel theory was given in a short 
note by one of the authors ^T|. However a systematic 
theory of classical Coulombic phase transitions including 
a comparison with numerical and experimental data was 
given only in the pioneering work of Michael Fisher and 
his coworkers |22j . These authors have also made an ex- 
tensive investigation of the state of art in this field p3| . 
Summarizing these results for classical systems we may 
say that theory and experiment (Monte Carlo data as 
well as measurements on electrolytes) are in quite good 
agreement. There is no doubt any more that a Coulom- 
bic phase transition in classical systems exists. We will 
not repeat here all the arguments but refer to the careful 
investigations of Fisher et al. (2^j2^] . 

The classical Coulombic transition is due to a balance 
between a hard-core repulsion and a Coulombic attrac- 
tion. We will show here that the PPT is a balance 
between quantum repulsion between point charges and 
Coulombic attractions. In this respect PPT is a kind of 
quantum variant of the classical Coulombic transition. 
If on one hand the existence of a Coulombic transition 
in ionic systems is now well confirmed, there is on the 
other hand no proof yet for the existence of a PPT in 
plasmas. At present the problem of the existence of a 
PPT is still open. From the point of view of the theory 
this is connected with the difficulties to derive an ac- 
curate equation of state for nonideal quantum plasmas. 
From the experimental point of view the difficulties are 
connected with the very high pressures where the PPT 
could occur. However as we will point out here there 
are now several experimental and numerical data which 
point to the existence of a PPT in real plasmas. 



II. DEBYE HUCKEL APPROXIMATION FOR 
CLASSICAL AND QUANTUM SYSTEMS 



tions in the regions where - according to the first es- 
timates - the phase transition is to be expected. As 
we have learned from the theory of phase transitions in 
neutral gases, already simple expressions as the van der 
Waals equation might give qualitatively correct results. 
We will show here that the Debye-Hiickel expressions for 
the thermodynamic functions together with the mass ac- 
tion law provide an approximation which might serve a 
a zeroth approximation for the description of Coulom- 
bic phase transitions. In their first work Norman and 
Starostin Q proposed for a first estimate of the criti- 
cal temperature, to use Debye-Hiickel-type expressions 
(1 — cX/td) for the chemical potential. Here A is electron 
thermal wavelength, 775 is a Debye radius, c ~ 0.1 is a 
numerical coefficient. The Debye-Hiickel term describes 
an effective Coulomb attraction, the (A/r£>)-term rep- 
resents in this approach an effective quantum repulsion 
due to the uncertainty principle. Norman and Starostin 
H used also the original expression of Debye- Hiickel for 
the chemical potential (l + cX/ru)^ 1 which has the same 
accuracy if (A/ro) tends to zero. 

We will give now a systematic discussion of the Debye- 
Hiickel approximation for the classical and quantum 
cases. We consider a binary Coulomb system with n+ 
positive ions (cations) and n_ negative charges (anions or 
electrons) per cubic centimeter n+ = n_. The density of 
neutrals is no- The total density is n = n + +no = ?i_+no. 
In the following we will use the plasma notations, i.e. 
we call the negative charges "electrons" and the posi- 
tive charges simply " ions" . In the Debye-Hiickel theory 
for charged hard spheres with the diameter a the excess 
chemical potential of the charges of species i is given by 
the simple expression 
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(1) 



The densities of free and bound particles are connected 
by the mass action law (Saha equation) 
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where K(T) is the mass action constant which is given 
by Bjerrum's formula or certain modifications J5|,Q in 
the classical case. The osmotic pressure is given by the 
formula 



P/kT = m 
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Here the Debye-Hiickel function <fi(x) is defined by 
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In order to describe Coulombic phase transitions one 
needs good approximations for the thermodynamic func- 



The easiest way to transfer the Debye-Hiickel approxi- 
mation to quantum plasmas is the replacement of the 
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hard sphere diameter by an effective quantum diameter 
of charges. According to the Heisenberg principle a point 
particle with the thermal momentum 



P 



8mkT 



corresponds to a wave packet of size 
P 4V2nmkT' 



Sx 



(5) 
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Identifying this with the effective diameter of the charges 
we get 



»m = I = 



8y/2nm e kT 



In this so-called Lambda approximation |7|jL0|| quantum 
effects are expressed by just one characteristic length, the 
thermal de Broglie wavelength A in combination with 
Debye-Huckel type approximations. The Lambda ap- 
proximation agrees very good with the exact quantum- 
statistical results for small densities for T < 10 5 AT |7|,[|. 
This approximation requires that the plasma is nondc- 
generate i.e. n*A 3 << 1 (where n* denotes the density 
of free particles). Furthermore only a small fraction of 
the atoms should be bound in molecules 02 « 0.5. The 
standard choice of the mass action constant for quantum 
plasmas is the Brillouin-Planck-Larkin expression 



K(T) = A 3 a{T) 



(8) 



where 



In the case of classical ionic solutions the parameter a is 
to be identified with the diameter of ions. In the case of 
dense quantum plasmas we have to introduce the quan- 
tum effective diameter a(T) = A/8. This way we have 
shown that the PPT is the analogue of the known clas- 
sical Coulombic transition in ionic solutions. At least in 
the simplified theory given above the classical transition 
in ionic solutions and the quantum PPT in dense plas- 
mas are in full analogy. In both cases the spinodal curve 
is given by 
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(7) where b = e /ksTa(T) and [i = e^/^T. In the quan- 
tum case we find the corresponding estimate for the crit- 
ical temperature 
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(14) 



In the case of electrolytic solutions or semiconductors the 
critical temperature is much lower due to the appearance 
of the dielectric constant in the denominator. Another 
effect which decreases T cr is the finite size of ions. In or- 
der to describe the size of the ion cores in alkali or noble 
gas plasmas we introduce an effective ion radius R, e.g. 
we have R ~ 1.69 A for Cs— ions. Further we replace 
A/8 by 



a(T)= A/8 + R 



(15) 



Assuming that R is temperature-independent we find 
now from Eq.(15) 



a(T) =J2 s2 ( e M~PE s ) - 1 + 0E.) 



(9) 



In this way the classical and the quantum case might be 
treated in the present rough approximation by the same 
procedure. The easiest way to check for stability of the 
system is the investigation of the sign of the derivative of 
the pressure. The region where 
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corresponds to thermodynamic instability. Equivalent is 
the condition 



drij 
dn 
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(11) 



which is valid only if atoms are treated as ideal gas. In- 
stability is restricted to temperatures below the critical 
one which is defined by 



T < T r/ , 
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T rr = 12570^ 



1 + (1 + {6R/7Ta B )) 1/2 



(16) 



where a B = h/e 2 m e is the Bohr radius. We see that a 
finite ion size may lead to an essential decrease of the crit- 
ical temperature to values far below lOOOOA'. For Cs— 
plasmas this estimate gives e.g. T cr ~ 6000^. Since 
this value is still far above the observed critical point 
T cr ~ 2000-RT we have to search now for additional ef- 
fects, which might reduce the critical temperature. 



III. COMBINATION OF VAN DER WAALS AND 
DEBYE-HUCKEL APPROXIMATION 

In the Debye-Huckel type approximation derived in the 
previous section the interaction of charges with neutrals 
and the neutral-neutral interaction was not taken into ac- 
count. This model might fail if the interactions with neu- 
trals have an essential influence. In the following we con- 
sider for concreteness a plasma system. We treated the 
neutrals as ideal particles in the simple approach given 
above which is analytically solvable. The system of mass 
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action laws might become rather difficult if one includes 
the interactions with neutrals or/and chemical equilibria 
between neutrals as e.g. between atoms and molecules. 
A better starting point for such more complicated cases 
is as a rule the expression for the free energy, including 
a minimization procedure. 

Combining the Debye-Hiickel model with a van der 
Waals expression and a contribution which takes into ac- 
count the polarization of neutrals we get for the free en- 
ergy 



F = n k B T 
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Here the last three terms denote the contributions from 
the short range van der Waals interactions and the contri- 
butions from polarization terms. The polarization terms 
are of particular importance for the description of alkali 
and mercury plasmas p5| , p6|] . 

Three free constants characterize now our model 

- the repulsion of atoms : B 

- the attraction of atoms : A 

- the strength of atomic polarizability : W 

In order to reduce the number of free parameters we 
assume in the following W = 0. We note that the free en- 
ergy in this form defines a model plasma with properties 
intermediate between van der Waals gas and a Debye- 
Hiickel plasma. The equilibrium composition is given by 
the nonideal Saha equation 



1 — Q 
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Here 



AI(a) 



e 2 Ky/a (1 — a)nB 

' k B T(l + Ka(T)Va) ~~ 1 - (1 - a)nB 
+ log(l - (1 - a)nB) 

+ (1 - a)nA - Wn(2a - 1) 



(19) 



is the lowering of the ionization energy. 

As well known, real gases with attractive interactions 
show a first-order phase transition described in the p — T 
plane by a critical point G\ and a coexistence line ending 
in C\. Assuming that the Coulombic terms are omitted 
in our expressions for the free energy a simple van der 
Waals type expression is obtained with the critical point 



TvdW ~ TikVb^ - Yb 



(20) 



The critical temperature is about 10 1 < T cr < 10 3 for re- 
alistic values of the parameters. For example T cr ~ 33 AT 



for hydrogen. The degree of ionization is practically zero 
in this region, i.e. there is no interference with the PPT. 
On the other hand, at least for hydrogen, the PPT oc- 
curs in a region where the number density of neutrals is 
rather low. We may assume therefore that the assump- 
tion a = 0, b = is justified in this case. Then, as shown 
above, the PPT may be treated in an analytical way In 
the Lambda-approximation we got in the previous sec- 
tion an analytical expressions for the critical point of the 
PPT 0,|l^]. By introducing the Bohr radius a B we may 
transform the expression for the critical temperature to 
the form 



Tppt ~ 8k B a B 



(21) 



Further the critical density of the free electrons may be 
written as 



(22) 



What happens in a system where both types of in- 
teractions are present? Evidently now besides the 
classical first-order phase transition typical for neu- 
tral gases a second first-order phase transitions due 
to Coulomb forces may appear. The second one 
which might occur only at rather high temperatures 
T > 10 4 A is the PPT. In Fig. 1 we give the pres- 
sure isotherm of hydrogen plasmas at T = 10000 A. 
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FIG. 1. The isotherm T = lOOOdft: for hydrogen plas- 
mas showing the PPT instability. 

In the unstable region of the isotherm the plasma is 
separated into two phases. The dense phase is highly 
ionized and the Coulomb interaction dominates over the 
van der Waals forces. The less dense phase is only weakly 
ionized and van der Waals interactions dominate. A 
more realistic calculation of the coexistence line for hy- 
drogen plasmas was given recently pfij ]. A coexistence 
pressure of p ~ llbGPa was found in the cited work for 
T = 10000A; the corresponding mass densities of the co- 
existing phases are p\ ~ 0.62(;/cto 3 and p2 — 0.82g/cm 3 . 

In order to investigate the qualitative influence of the 
van der Waals forces we studied a model plasma with 
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increasing values of the parameters a and b. Fig. 2 
shows the isotherm T — 0.07 a.u. and a = b = 30 a.u.. 
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FIG. 2. The isotherm T = 0.07 a.u. for an alkali-like 
model plasma with one transition of mixed PPT-van 
der Waals type (van der Waals parametersof the model: 
a — 30 a.u., b = 30 a.u.). 

In spite of the fact that these values are rather large we 
cannot detect a separate van der Waals wiggle. Both at- 
tracting forces support each other and lead to one phase 
transition approximately in the same density region as 
the PPT. This is the situation we observe in nature for 
alkali plasmas. At further increase of the values of the 
van der Waals parameters to a = b = 100 a.u. we ob- 
serve two separate van der Waals loops as shown in Fig. 
3. This is the situation we observe for real noble gas 
plasmas. We are to emphasize, that Fig. 2 and Fig. 3 
have just illustrative character; in reality we do not see 
two wiggles in one isotherm due to the drastic difference 
between the critical temperatures. 

Summarizing these findings we may state that in de- 
pendence on the values of the van der Waals parameters 
we may obtain either a phase diagram with two first order 
transitions (hydrogen and noble gas plasmas) or a phase 
diagram where both transitions fuse to just one (alkali 
plasmas). The existence of a phase diagram including 
a van der Waals type transition and a separate metal- 
insulator phase transition was discussed for the first time 
in 1944 by Landau and Zeldovich Jjj]. We repeat that the 
existence of two separate phase transitions requires the 
validity of the inequalities 



T v dw « T ppt ;p v dw « Pppt 



(23) 



This is fulfilled for hydrogen and for noble gas plasmas. 
In the case of alkali plasmas the van der Waals and 
the PPT transition have about the same critical den- 
sity and temperature. Therefore both transitions fuse 
just to one with very specific properties HJ2^M|B7| . 
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FIG. 3. The isotherm T = 0.07 a.u. for a noble 
gas-like model plasma with a PPT and a separate van 
der Waals transition (van der Waals parameters of the 
model: a — 100 a.u., b = 100 a.u.). 



IV. DISCUSSION OF PLASMA PHASE 
TRANSITIONS 



Considerable efforts were applied to obtain better 
plasma thermodynamic functions for larger range of the 
nonideality parameter. We will not go into the details of 
the theory here which are explained elsewhere |7]|],[l(||llJ . 
The critical temperatures derived from more refined ver- 
sions of the theory lay in the region T cr = 14000- 19000^ 
[fflp||Tl]] . There is a scatter in the values of T cr , p cr and 
n cr , since the procedure of getting the thermodynamic 
functions for large nonideality is not a unique one. In- 
stead of going into details on the several theoretical ap- 
proaches we concentrate here on the physics and the re- 
lation to experiments, to numerical simulations and to 
related phenomena. 

Let us start with the experimental aspects. The exis- 
tence of a PPT in dense plasmas is still less clear than 
in classical Coulombic systems. This is mainly due to 
the fact that pressures predicted for the coexistence line 
are above IMbar what is still hardly reachable in ex- 
perimental situations. However in the last time more 
and more experiments were performed which cover the 
region of interest |9||^|^,||,||||,|||| . Experiments 
with shock-compressed hydrogen and deuterium plasmas 
have verified that around 140GPa and 3000-ftT a transi- 
tion to a highly conducting state occurs p^ , Bl| . This is 
approximately the pressure region where the coexistence 
line of the PPT is expected |^8) . A corresponding behav- 
ior of conductivity data has been reported 30 1 . Pressure 
dissociation and ionization, which almost do not depend 
on the temperature, become a dominant factor at 10 4 JsT 
and below in hydrogen and deuterium fluid studied ex- 
perimentally at high pressures |3l|,p4 35 39]. 



5 



In order to compare theory and experiment in a more 
quantitative way refined treatments of the phase tran- 
sition induced by pressure dissociation and ionization 
were given RMO,|l3|,[l4 28 1 . Most of the possible species 
and various interactions were taken into account to cal- 
culate the equations of state. There is a certain dis- 
cussion whether the model equations of state reproduce 
reasonably well the recent experiments. However the 
theory predicts in all variants a phase transition, with 
T cr = 14000- 19000.fr . Excited states do not contribute 
significantly to the hydrogen partition function in the in- 
termediate range of temperatures. 

Path integral Monte Carlo calculations and wave 
packet dynamics also give some hints to the existence of 

the plasma phase transition QHIllllimSH ' F ° r 
example the validity of the thermodynamic expressions 
used was tested by comparison with molecular simulation 
and experimental data for nonideal plasmas and limiting 
expansions for weakly nonideal plasmas fjjjfl . The prob- 
lem continues to attract remarkable attention, for exam- 
ple at the recent (June 2-7, 2002) "International Confer- 
ence on Warm Dense Matter and FEL Experiments" in 
Hamburg, Filinov et al [ ^4||45| ] claim that they "present 
results of direct path integral Monte Carlo simulations 
which, for the first time, provide first-principle support of 
such a phase transition in dense hydrogen" . It should be 
noted, however, that these particular PIMC simulations 
yield rather low energy values for the "droplets" that 
appear in the simulations. Kohanoff and Scandolo |4(| 
discussed on this workshop the same problem in a softer 
manner and presented their results of ab initio molecular 
dynamics simulations. We mention these works as part 
of the ongoing effort to simulate the PPT. 

Let us discuss now in brief the phase transitions ob- 
served in alkali-plasmas. As pointed out in previous 
sections the theory predicts for this case that the PPT 
and the van der Waals transition fuse just to one transi- 
tion. Near to the critical point the degree of ionization is 
low but different from zero. Coulombic interactions are 
present and influence the critical behavior, in spite of the 
fact that they do not play the dominant role. The area 
of parameters close to the critical point of cesium was 
carefully investigated [^5| . It was shown that near to the 
critical point the ionization landscape changes quickly 
and has a quite complicated structure. The theoretical 
predictions which take into account the charge-atom in- 
teraction (polarization effects) are in rather good agree- 
ment with the data 26 14|. 

The existing experiments for argon plasmas based on 
shock wave data are also in satisfactory agreement with 
the theory ggg. 

We will discuss now the problem of critical behavior 
and the critical indices. For the classical case much atten- 
tion has been devoted to this problem 
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For the 

PPT where quantum effects influence the critical point, 
the problem of critical indices is still open. So far no 



experiments are available which give reliable information 
about the critical properties of hydrogen and noble gas 
plasmas. Experimental studies of classical electrolyte so- 
lutions revealed that Coulomb liquids characteristics dif- 
fer from those of simple liquids. Even if crossover was ob- 
served from classical scaling laws to Ising scaling, it takes 
place much closer to the critical point. The Ginzburg pa- 
rameter, which characterizes the size of the Ising region, 
turns out to be of 1-2 and more orders of magnitude less 
than for simple liquids. So there is a remarkable dif- 
ference in crossover range between simple and Coulomb 
liquids |23|,|^|[50| . Experimental study of crossover phe- 
nomena in real plasmas and in particular in cesium plas- 
mas might help to find out if the critical point is of the 
plasma or gas-liquid kind. Additional factors, which do 
not exist in electrolytes, namely, quantum effects and 
charge-atom interaction can influence the crossover phe- 
nomena in PPT. 

Michael Fisher (49) formulated clearly the challenge 
to the critical phenomena theory: "Critical behavior is 
thoroughly understood in Ising models alias lattice gases: 
But how far does that help our insight into the critical 
region of continuum models and real fluids? In particu- 
lar, what causes some ionic solutions to exhibit van der 
Waals or "classical" critical exponents, while other, so- 
called " solvophobic" systems are of Ising - type?". It 
is evident that the difference between properties of var- 
ious liquids under similar external conditions can be at- 
tributed principally only to the interaction potentials be- 
tween the particles. However the Kadanoff- Wilson theory 
of critical phenomena, which is based on the universal- 
ity hypothesis, neglects interaction potential peculiari- 
ties and appeals only to large-scale properties, such as 
system dimension and Hamiltonian symmetry. So other 
approaches are needed. 

First of all Fisher points to the integral equation ap- 
proach. But most of the standard approaches to the 
theory of critical phenomena based on integral equa- 
tion method fail (4^] . Though both the Kadanoff- Wilson 
theory and the integral equation approach are on equal 
rights, the first is based on global treatment whereas 
the latter is related to the local Ornstein-Zernike equa- 
tion. For this reason Fisher puts the integral equation 
approach on the first place in his hope to solve the prob- 
lem of Coulomb criticality. However there is a specific 
difficulty in the local approach which is connected with 
the choice of closing relation in the integral equation ap- 
proach. This difficulty was attacked recently by Mar- 
tynov fl|. He suggested an integral equation which de- 
pends on the interaction potential and permits in prin- 
ciple to define the critical indices and their dependence 
on the interaction potential. Martynov [ [si] applied his 
local approach to the system of charged particles. The 
critical indices obtained differ almost as much as twice 
from van der Waals and Kadanoff- Wilson predictions. 
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However the results cannot be applied directly to elec- 
trolyte solutions since they do not take into account sol- 
vent molecule influence, neither it can be applied directly 
to dense plasma, since quantum effects are not taken into 
account. 

Let us finally discuss several special problems con- 
nected with Coulombic phase transitions. At first we 
mention the magnetic field influence. It was shown re- 
cently that strong magnetic fields increase the critical 
temperature |52| , |l2| ] . Another important effect is con- 
nected with nonequilibrium phenomena. It was noted 
p3[ that nonideal plasmas are usually generated in non- 
equilibrium state with respect to plasma wave excitation. 
This might complicate the theory Jl4| and the analysis 
of experiments. 

At third we would like to discuss a new class of phe- 
nomena observed in Cs-plasmas Though the stan- 
dard cesium phase diagram is well known |2^| , Recently 
a new phase state of Cs was observed by Holmlid et 
al [5^j55| , |5^| . In experiments with thermoionic convert- 
ers clusters were observed first at 1300if. These cluster 
were cooled and associated in microdroplets. The micro- 
droplets formed were fixed and measured at 7QK. The 
number density obtained was about 1018 atoms/cm 3 . 
Holmlid |54|] considered his results as an observation of 
the Rydberg matter predicted by Manykin et al p5| , p6[ . 
Generally speaking Rydberg matter exists only at OK 
temperature. At higher temperatures Manykin's ap- 
proach can and should be combined with our non-ideal 
plasma treatment. In fact both approaches enhance each 
other. The phenomenon might be connected with the 
idea of isolated segment of metastable nonideal plasmas 
which was introduced into the plasma thermodynamics 
long ago |57| . Supercooled nonideal plasma states were 
a direct consequence of the PPT [Q. Manykin's Ry- 
dberg matter f35|j56f , Holmlid's microdroplets and 
plasma phase transitions OME^ are various facets of 
Coulombic phenomena. The existence of isolated regions 
of metastable nonideal plasmas is not an exceptional fea- 
ture of cesium. Several authors |57],|5^,|55|,^6| suggested 
that this might be a more general feature. In fact Holm- 
lid observed his clusters in many other substances as 
well. The isolated region of metastable nonideal plasmas 
is possibly an inherent part of the phase diagram. This 
part complements the standard diagram and is superim- 
posed on it. 



the PPT is an analogue of the known Coulombic transi- 
tion in classical ionic solutions. In both cases the mecha- 
nisms are the same, the transition is due to the interfer- 
ence between a Coulombic attraction of the charges and 
repulsive forces. The classical repulsion of ions at small 
distances is in the quantum case replaced by Heisenberg 
and Pauli repulsion effects. We have shown above for 
dense gaseous plasmas that in dependence on the values 
of the van der Waals parameters of the neutral compo- 
nent we may obtain either a phase diagram wi two first 
order transitions (hydrogen and noble gas plasmas) or a 
phase diagram where both transitions fuse to just one (al- 
kali plasmas). In spite of considerable efforts there is no 
final proof yet that the PPT really exists, however new 
theoretical and experimental results seem to confirm this 
hypothesis. One strong theoretical argument in favor of 
the PPT is the analogy to the Coulombic phase transition 
in classical systems. For ionic solutions theory and ex- 
periment (Monte Carlo data as well as measurements on 
electrolytes) show a corresponding Coulomb transition. 
Therefore there is no doubt any more that a Coulombic 
phase transition in classical systems exists [22 2^ , (49| . For 
quantum systems as e.g. hydrogen plasmas the problem 
of the existence of a PPT is still open. The main reason 
for still unsolved problems are 

• the theory of the PPT requires an accurate theory 
of dense quantum plasmas, and 

• the pressures where the PPT occurs are in the 
Mbar region what leads to very difficult experimen- 
tal problems. 

We discussed here several results of theory and experi- 
ment which point to the existence of a PPT in real plas- 
mas and discussed the conditions and parameter regions 
where the PPT is expected. 
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V. CONCLUSIONS 

This work is devoted to plasma phase transitions 
(PPT) which are determined by attractive Coulombic 
interactions and quantum repulsion. The PPT occurs 
in the region of high pressures and temperatures, where 
matter is at least partially ionized. At least in principal 
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